Spectral functions for strongly correlated 5f-electrons 
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We calculate the spectral functions of model systems describing 5/-compounds adopting Cluster 
Perturbation Theory. The method allows for an accurate treatment of the short-range correlations. 
The calculated excitation spectra exhibit coherent 5/ bands coexisting with features associated 
with local intra-atomic transitions. The findings provide a microscopic basis for partial localization. 
Results are presented for linear chains. 



I. INTRODUCTION 

Electronic correlations are strongly evident in actinide 
intermetallic compounds. The highly complex phase dia- 
grams with novel sometimes enigmatic ordered states re- 
flect the sensitivity to small changes in external control 

axameters like temperature, pressure, magnetic fields 
. An undispensible prerequisite for an explanation of 
the observed anomalies is a microscopic understanding 
of the strongly interacting 5/ electrons which partially 
preserve atomic-like character. Occupying the partially 
filled 5/ shells according to Hund's rules leads to mag- 
netic moments. The lifting of rotational symmetry in a 
crystal by the crystalline electric field (CEF) and the hy- 
bridization with the delocalized conduction states of the 
outer shell electrons leads to a large number of low-energy 
excitations. 

Of particular interest are the heavy fermion phases 
where the low-energy excitations correspond to heavy 
quasi-particles. Their enhanced effective masses m* are 
reflected in enhanced values of the Sommerfeld coeffi- 
cient and the Pauli spin susceptibility. The heavy quasi- 
particles have been observed by de Haas-van Alphen (dH- 
vA) experiments in a number of compounds. The experi- 
ments unambiguously conflrm that some of the U 5/ elec- 
trons must have itinerant character. It has been known 
for quite some time that the 5/-states in actinide in- 
termetallic compounds cannot be considered as ordinary 
band states. Standard band structure calculations based 
on the Local Density Approximation (LDA) fail to repro- 
duce the narrow quasiparticle bands. On the other hand, 
the predicted bandwidths are too small to explain pho- 
toemission data 0,01 ■ These shortcomings are a direct 
consequence of the inadequate treatment of local corre- 
lations within ordinary electron structure calculation. 

Increasing experimental evidence points towards a two- 
fluid model implying the co-existence of both localized 
atomic-like and itinerant band-like 5/ states. For U 
heavy fermion compounds the photoemission spectra 
usually display a two-peaked structure in the 5/ emis- 
sion ^ ^] . The peak at 1 eV binding energy results from 
the localized 5/ states while the narrow peak at the Fermi 
level is attributed to 5/ derived itinerant quasi-particles 
0. Concerning the low-energy excitations it has been 
shown recently that the dual model allows for a quan- 
titative description of the renormalized quasi-particles - 



the heavy fermions - in UPd2Al3. The measured dH-vA 
frequencies for the heavy quasiparticle portions as well 
as the large anisotropic effective masses can be explained 
very well by treating two of the 5/ electrons as local- 
ized Finally, the co-existence of 5/-derived quasi- 
particles and local magnetic excitations in this com- 
pounds have been conflrmed by recent neutron scattering 
experiments Q. Theoretical studies aiming at an expla- 
nation of the complex low-temperature structures lay em- 
phasis on the partitioning of the electronic density into 
localized and delocalized parts Considering the 
success of the dual model in describing the low-energy 
behavior of actinide heavy fermion compounds we have 
to explain how the splitting of the 5/ shells arises and 
propose means of identifying the underlying microscopic 
mechanism. 

The dual model as sketched above is an effective Hamil- 
tonian for the low-energy regime. In the spirit of Wilson's 
renormalization group it should be obtained by integrat- 
ing out processes at higher energies. The central the- 
oretical task is to identify the microscopic mechanisms 
which cause the 5/ band widths to renormalize to al- 
most zero for certain orbital symmetries while staying 
flnite for others. In the present paper we concentrate on 
the role of intra-atomic correlations, the motivation be- 
ing as follows: Recent model studies for the ground state 
properties of small clusters show that intra-atomic cor- 
relations as described by Hund's rules may strongly en- 
hance anisotropics in the kinetic energy and thus lead to 
an orbital-selective Mott transition in 5/-systems 
This scenario closely parallels the one postulated re- 
cently to explain the behavior of transition metal oxides 
[l3,[lJ,llEllfl|- Here we discuss how intra-atomic corre- 
lations affect the single-particle spectral functions. The 
conjecture is that the above mentioned enhancement of 
anisotropics in the kinetic energy can explain the dual 
nature of the 5/ electrons. In a real material this en- 
hancement would imply anisotropics far beyond these 
predicted by a standard electronic structure calculation. 

The concept of correlation-driven partial localization 
in U compounds has been challenged by various authors 
(see e. g. ^3|)- The conclusions are drawn from the fact 
that conventional band structure calculations within the 
Local Density Approximation (LDA) which treat all 5/- 
states as itinerant can reproduce ground state properties 
like Fermi surface topologies, densities. The calculation 
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of ground state properties, however, cannot provide con- 
clusive evidence for the delocalized or localized charac- 
ter of the 5/-states in actinides. First, the presence of 
localized states can be simulated in standard band cal- 
culations by filled bands lying (sufficiently far) below the 
Fermi level. Second, the Fermi surface is mainly deter- 
mined by the number of particles in partially filled bands 
and the geometry of the lattice which affects the disper- 
sion of the conduction electrons. A change in the number 
of band electrons by an even amount does not necessarily 
affect the Fermi surface since a change by an even num- 
ber may correspond to adding or removing a filled band. 
As such, the Fermi surface is not a sensitive test of the 
microscopic character of the states involved. 

Following up this track of thought we would like to 
present here a qualitative discussion of strongly corre- 
lated 5/ systems emphasizing properties of the spectral 
functions which would allow an experimental discrimina- 
tion between correlation-driven partial localization and 
a usual filled-band scenario. We lay emphasis on the 
role of orbital degeneracies and on intra-atomic correla- 
tions which compete with the (anisotropic) kinetic en- 
ergy. The features chosen for discussion are the distribu- 
tion of spectral weight, the formation of coherent bands 
and the opening of an excitation gap. 

We begin the discussion in Section ^ with a brief ex- 
position of the model Hamiltonian and the techniques 
emphasizing the general structures involved. The numer- 
ical procedure and the parameters adopted in the actual 
calculation are the presented in Section UTTI In the sub- 
sequent Sections IIVI and EI we delineate details of the 
models and the results obtained. The conclusions are 
summarized in Section IvTI 



II. MODEL HAMILTONIAN AND 
COMPUTATIONAL METHOD 



For the adoption of a model describing the 5/ elec- 
trons in actinide based heavy fermion systems we rely 
on experience with ab initio electronic structure calcu- 
lations. In these materials, the direct overlap of the 5/ 
wave functions at neighboring sites is anticipated to be 
rather small due to the large U-U distance. Detailed 
studies suggest that the 5/ states acquire their dis- 
persion by hybridization with high-lying empty conduc- 
tion states. We model these processes by introducing 
weak effective transfer integrals. The interplay between 
intra-atomic Coulomb interaction and anisotropic kinetic 
energy in 5/-systems is described by the simple model 
Hamiltonian 



H = H, 



band ^ -'-'Coul 



Ha 



(1) 



^Sli (°)^J-3 (")^J"^4 («)• (2) 

is written in terms of the usual fermionic operators (a) 
(cj^(a)) which create (annihilate) an electron at site a 
in the 5/-state with total angular momentum j and 
z-projection jz- Considering the fact that the spin- 
orbit splitting is large we neglect contributions from 
the excited spin-orbit multiplet j = 7/2 and adopt the 
j — j'— coupling scheme. The Coulomb matrix elements 
U. 
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are given in terms of the usual Clebsch-Gordan coeffi- 
cients C \ and the Coulomb parameters U j. 

The kinetic energy operator in the band term i?b„,d 
describes the hopping between all pairs neighboring sites 

(ah) 

• (4) 

We assume the transfer integrals tj^ to be diagonal in 
the orbital index jz- This simplification is justified for 
a 1-dimensional system since the cylindrical symmetry 
allows one to select a common quantization axis for all 
sites. We would like to stress that although the hopping 
term is assumed to be diagonal in the orbital indicies the 
corresponing channels are coupled through the Coulomb 
term. Finally, we account for the orbital energy e/ which 
determines the /—valence of the ground state. 
The single-particle spectral function 
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with 



The local Coulomb repulsion 



yields the probability for adding {us > 0) or removing 
(w < 0) an electron with energy w in a state characterized 
by momentum k and angular momentum projection jz 
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Figure 1: Representation of a one-dimensional chain as a 
super-lattice of identical L-site cluster. The orbital-dependent 
transfer integrals between neighboring sites within a cluster 
and between adjacent clusters are t = (tj^^_c^/2, ■ ■ ■ 



and to = (t, 



-5/2, 



respectively. 



=5/2j 



% with energy 



{N) 



to the iV-particle ground state 
The states with iV ± 1 and their energies are denoted by 
^{N±i}\^ and En^^^\ respectively, and L is the number 
of sites. 

The importance of the spectral functions stems, first, 
from the fact that they can be observed in photoemission 
and inverse photoemission experiments and second, that 
the orbital-projected expectation value of the kinetic en- 
ergy per site 



/ + 00 1 
dLof{Lo)-J2e,Ak)A,Ak,co) 
-oo ^ 1, 



(7) 



as well as the orbital-projected momentum distribution 
function 



rij^ (k) 



-\-oc 



duj f{uj)AjAk,uj) 



(8) 



and the orbital-projected Density of States (DOS) 

k 

can be expressed through them. The dispersion of the 
band energies ej^ (fc) is given by the kinetic energy Eq. 
Q , L denotes the number of lattice (cluster) sites while 
/ is the usual Fermi function. 

We determine the spectral function Aj__{k,uj) from the 
single particle Green's function Gj^{k,uj), 



AjAk,(^) ^ lim Im Gj,(fc, w -f i77) (10) 

which, in turn, are calculated applyi ng C luster Pertur- 
bation Theory (CPT) [UliS HI El 0] ■ The method 
which has been successfully appHed to a wide variety of 
many-particle models proceeds by dividing the infinite 
periodic lattice into identical finite clusters as illustrated 
in Figure ni 

The Hamiltonian Q takes the form 



and n. The Green's functions of the finite clusters 



(Gf)(.)) 



ab 



(12) 



are calculated exactly thereby fully accounting for the 
complex dynamics of the strongly correlated system. 
Here we use the locator representation where the indices 
refer to the sites in the cluster. In a second step, the 
inter-cluster hopping VJ^j^ , is treated as weak perturba- 
tion which yields the Green's function of the super-lattice 
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The variation with the wave vector Q which is taken from 
the first Brillouin zone of the super-lattice is introduced 
by the inter-cluster hopping matrix Vj^{Q). Its matrix 
elements are obtained from V^^^ by performing a Fourier 
transformation with respect to the cluster indices m and 
n. Following the translational invariance of the orig- 
inal lattice is explicitly restored yielding the Green's func- 
tion 



(14) 



Note that the CPT is exact in three limiting cases: van- 
ishing hopping terms, in an uncorrelated system and for 
infinite cluster size. The approximation is controlled 
most effectively by considering different cluster sizes L 

Before we turn to a detailed description of the calcula- 
tions we would like to mention the following caveat. The 
formulation of CPT as summarized above implicitly re- 
lies on the assumption that the ground states of the inter- 
acting electron system are non-degenerate. This propo- 
sition, however, may not be satisfied for the building 
blocks of the super-lattice. The problem arises trivially 
for clusters containing an odd number of electrons where 
all states are at least two-fold degenerate. Of course one 
could try to circumvent the problem by choosing appro- 
priate cluster sizes or by varying the band filling, i. e., 
the number of particles. This strategy which may work 
for simple systems does not provide a solution in orbitally 
degenerate systems. Due to the local orbital degrees of 
freedom the complexity increases rapidly with the num- 
ber of sites or particles imposing serious restrictions on 
models. As we attempt at a description of homogeneous 
phases we separately calculate the Green's functions for 
the states of the ground state manifold and subsequently 
average the spectral function. 



^ = E^™+ E KT," cL,.c„6.., (11) 

m m,n,a,b,jz 

where i?^ refers to cluster m while V^Jl describes the 
hopping between sites a and b on the adjacent clusters m 



III. COMPUTATIONAL DETAILS 

We begin by specifying the model parameters entering 
the Hamiltonian Eq. 1^ . The high-energy scale is set by 
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the Coulomb repulsion between two 5/ electrons at the 
same site. The parameters Uj are chosen according to 
the following considerations. First, it is well known, that 
the isotropically averaged Coulomb repulsion in a metal 
is strongly reduced as compared to its value in an atom. 
The reduction is a direct consequence of screening by 
the itinerant conduction electrons which, in turn, implies 
that the actual value depends upon the chemical environ- 
ment of the correlated ion under consideration. We do 
not attempt at an ab initio calculation this quantity but 
rather leave it as a parameter, Uat, whose value can be 
estimated from the positions of the valence peaks in pho- 
toemission and inverse photoemission. The intra-atomic 
correlation which are the focus of interest in the present 
paper involve the anisotropic parts of the Coulomb inter- 
action. The latter give rise to the multiplet structure. It 
is important to note that the corresponding interactions 
are (usually) not screened and hence retain their atomic 
values. As we expect the anisotropic Coulomb parame- 
ters to be rather robust with respect to changes in the 
chemical composition or to pressure we use fixed values 
for the differences 



AU4 = Uj=4 - Uj^o = -3.79eV 

AU2 = Uj=2 - Uj^o = -2.72eV . (15) 

which were determined from an ab initio calculation for 
UPts ^j. The value for the orbital energy £/ in Eq. 
Q is fixed by the requirement that low energy states are 
formed by linear combinations of the ionic and Z'^- 
configurations which are energetically (almost) degener- 
ate. The intermediate valent ground state with ri/ ~ 2.5 
allows for low-energy valence transitions and the forma- 
tion of bands also in the strong-coupling limit. 

Finally, the values of the transfer integrals tj^ are var- 
ied in the strong-coupling regime. As in Ref. [I2l we 
choose ti/2=t5/2 in the figures. 

We next turn to the calculation of the cluster Green's 
functions G^- (z) which account for the local correlated 
dynamics. We should like to mention that the use of the 
open boundary conditions required by the CPT method 
considerably increases the numerical effort over that en- 
countered in finite clusters with periodic boundary condi- 
tions. A further difficulty arises from the fact that the or- 
bital degrees of freedom may lead to a rather high density 
of low-energy states. Keeping in mind these complica- 
tions the Green's functions for clusters with two and three 
correlated sites are calculated using Lanczos continued- 
fraction technique To avoid convergence problems 
due to the near degeneracy of the lowest excited states in 
the 5/ systems we use the JDQR algorithm instead 
of the usual Lanczos method to obtain the ground state. 

For purposes of comparison we present alongside the 
CPT spectra the test cases of finite clusters using periodic 
boundary conditions. In this case, the calculation of the 



cluster Green's function Gj^ (z) simplifies to 

= i E e''^('^-'')G(f'(fc,,z) (16) 

where 

GfJiKz) = Gf:l{k,z) + Gf^l{k,z) 

G5fl(fc,^) = (*oli(fc)j^^^^— ^c,.(fc)|*o) (17) 

The discrete value kt = ;£ = 0,...,L-1 label 
the single-particle eigenstates of the L-site cluster. The 
Green's functions are evaluated for complex frequencies 
z — UJ + ir] with 77 > 0. The imaginary part leads to ad- 
ditional Lorentzian broadening of the spectral functions 
and may affect the behavior of the integrals like the ex- 
pectation value of the kinetic energy or the momentum 
distribution function. All these quantities are increas- 
ingly difficult to compute as ry decreases. For the actual 
calculations we choose values for rj of the order of 10~^ 
of the relevant hopping integrals. For the plots a rather 
large value of 77 = 0.03 is chosen to allow for a good 
visualization. 

The absolute values of Uj are very large compared to 
the hopping integral and thus we shall limit ourselves to 
the subspace of and configurations for the ground 
state which has shown to be a good approximation. This 
leads to a considerable reduction of Hilbert space dimen- 
sion and less computational costs. The actual size of the 
Hilbert spaces in the reduced systems for the two site 
cluster filled with five electrons are 600 states and 220 
respectively 400 for four and five electrons. For the clus- 
ters with three sites one obtains 18 000 dimensions for a 
filling with eight electrons and 8 000 respectively 13 500 
for seven and nine electrons. The matrix representations 
of the Hamiltonian are then stored as sparse matrices. 



IV. SPECTRAL FUNCTIONS OF FINITE 
CLUSTERS 

A. Atomic limit 

We begin by discussing the spectral function of our 
model in the atomic limit neglecting hopping between 
the sites. The results provide a first qualitative insight 
into the behavior in the strong-coupHng Hmit. Due to 
the rotational invariance of the Coulomb interaction Eq. 
(01 the spectral functions Aj^ (uj) will not depend on 
the magnetic quantum number jz- The lowest-energy 
states of an intermediate- valent system in the zero-band 
width limit are given by products of atomic /^- and Z'^- 
configurations whose energies are assumed to be (almost) 
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PES BIS 

Position / eV -1.06 0.00 0.00 1.60 2.73 
Spectral weight 0.054 0.196 0.218 0.032 0.083 

Table I: Positions and spectral weight of the atomic transi- 
tions in uranium from numerical calculations. For < eV 
the contributions are from the hole propagation and for cu > 
eV from the particle propagation. 

degenerate. The corresponding spectral functions are ob- 
tained in close analogy to the classical work by Hubbard 
[2^ . The zero-configuration width approximation which 
neglects intra-atomic correlations leads to a characteristic 
three-peak structure. The valence transitions 
and occur at large energies and, concomitantly, 

do not affect the low-temperature behavior. The latter 
is determined by the low-energy peak resulting from the 
transitions ^ within the /^- and /"^-configurations. 
This peak is a direct consequence of the intermediate- 
valent ground state. The strong correlations present im- 
ply a substantial transfer of spectral weight from the 
low-energy part to the high-energy regime. The weights 
of the peaks can be estimated from combinatorial con- 
siderations. The weight Z{f'^ f^) of the transition 
P p equals the probability that a state with a given 
jz is occupied in that p contribution of the mixed- valent 
ground-state. Following these lines one finds: 

Z{f Z{f -^P)^\ (18) 

The central focus of the present paper is the evolu- 
tion of the low-energy peak whose spectral weight sums 
up to 7/12. Intra-atomic correlations which are usu- 
ally described by Hund's rules further reduce the spec- 
tral weight of the low-energy excitations involving va- 
lence transitions , J — 9/2) ^ |/^,J = 4) between 
the ground state multiplets. The corresponding spec- 
tral weights can be expressed in a straightforward way 
in terms of the usual Clebsch-Gordan coefficients and 
the reduced matrix elements. In addition to the cen- 
tral peak one finds transitions to the excited multi- 
plets |/3,J-9/2) ^ |/^J = 2) and \P,J^4) ^ 
\P,J — 3/2) , \f^,J — 5/2) occurring at the corre- 
sponding multiplet excitations. The excitation energies 
as well as the corresponding spectral weights are listed 
in Table U We should like to mention that there is no 
transition into the excited multiplet | , J = O) . 

B. Results for finite clusters 

We next turn to the modifications introduced by the 
hopping between different sites which acts as a weak per- 
turbation in the strong-coupHng Hmit. Its principal ef- 
fect is to remove the rotational symmetry of the isolated 



atom and, consequently, remove the high degeneracy of 
the ground state. As the total angular momentum ceases 
to be a good quantum number the ground state may con- 
tain admixtures from excited J-multiplets. The model 
Hamiltonian as specified by Eqs. QJ-Q, however, con- 
serves the z-component of the angular momentum allow- 
ing us to classify the eigenstates with respect to Jz ■ The 
magnetic character of the ground states depends upon 
the transfer integrals as shown in Ref. These mod- 

ifications affect the spectral functions Aj^{ki,uj) of an 
L-site cluster where the discrete set of quantum num- 
bers kg = 0, . . . , {L—1)^ labels the single-particle eigen- 
states. First, the positions of the peaks corresponding to 
the transitions between different valence states may ac- 
quire dispersion. Second, spectral weight may be trans- 
ferred from the low-energy regime to high-energy satel- 
lites. Third, transitions to excited multiplets not present 
in the atomic limit may appear and, finally, the overall 
spectral weight may be redistributed among the differ- 
ent j,2-channels. These issues will be addressed in the 
following section. 

A first qualitative understanding can be gained by 
considering the spectral functions of a two-site cluster 
where simple approximate expressions are found for the 
ground states in limiting cases. Throughout the discus- 
sion we restrict ourselves to the strong-coupling regime 
\tj^\ <C |AJ7j| where the ground state is mainly built 
from products of local p configurations \p;4:Jz) and 
local p configurations |/'^;9/2J^). Finite hopping be- 
tween the two sites a and b splits the manifold in first or- 



j=±3/2 j=±1/2 j=±5/2 

'z 'z 'z 




-0.2 0.2 -0.2 0.2 -0.2 0.2 

M / eV CO / eV w / eV 



Figure 2: Variation with wave number of the orbital-projected 
spectral functions calculated for a two-site cluster with five 
electrons and the hopping parameters 1^/2 = t and ti/2 = 
^5/2 = i' ■ The full line and the dotted line refer to fe = 
and fc = TT, respectively. The spectra are obtained from aver- 
aging the Green's functions over the degenerate ground state 
manifold. Spectral weight is transferred to local excitations 
(valence transitions and transitions into excited atomic mul- 
tiplets) which are not displayed here. The Lorentzian broad- 
ening is = 0.03. 
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der perturbation theory. We consider the case of strongly 



anisotropic hopping |i3/2| 
and |<3/2| = |t| « 



\t\ » \t'\ = 1*5/2! 



1/2 1 

^5/2 1 = 1^1/2! as well as isotropic 
hopping t3/2 = t = t' — = ^1/2- The variation with 
t, t' of the spectral functions is displayed in Figure 121 
Let us first consider the behavior along the isotropic 



line t 



3/2 



t = t' ^ t 



5/2 



Cl/2- 



Due to the rotation- 



ally invariance, the magnitude of the total angular 
momentum J = J{a) + J(a) provides a good quan- 
tum number. It has been shown previously that the 
ground state is six-fold degenerate. It is a linear combina- 
tions of the J"^ eigenstates \{^A)JJ^) and \{^%)JJ2z) 
for J = 5/2 which can be expressed in terms of the 
above-mentioned product states \p] \ Jz — Jz) \P]'^Jz) 
and \p;AJz) \ f^; \ Jz — Jz)- The gain in kinetic energy 



1^1 



■ 4 J J, 



|t|f!|r(:^)l 



(19) 



given in terms of the reduced matrix element 



' yI and the 6j-symbol 



riJ) = (- 




(20) 



for = 5/2 results from the formation of symmetric 
(antisymmetric) combinations the above-mentioned 
eigenstates. The isotropy of the ground state is refiected 



j = ±3/2 



j = ±1/2 



= ±5/2 




-0.2 0.2 -0.2 0.2 
(o/eV co/eV 



Figure 3: Variation with wave number of the orbital-projected 
spectral functions calculated for a three-site cluster with 
eight electrons and the hopping parameters t-j/2 = t and 
ti/2 = t5/2 = t' ■ The full line and the dotted line refer 
to fe = and fe = respectively. The spectra are ob- 

tained from averaging the Green's functions over the degen- 
erate ground state manifold. Spectral weight is transferred 
to local excitations (valence transitions and transitions into 
excited atomic multiplets) which are not displayed here. The 
Lorentzian broadening is 77 = 0.03. 



t = 
f . -0.1 



t.-O.I 
f = 



t = 
f . -0.1 



(a) 



(b) 



Figure 4: Ratio Tj^/tj.^ of the two site cluster (left) and three 
site cluster (right) with tj^ < along the line tz/2 ~ t and 
ti/2 = ~ t' connecting the values below the figures. 



in the isotropic spectral functions which exhibit peaks at 
the energies 



±^ = ±0.077 eV 

14 1386 



(21) 



Due to the strong correlations we find two peaks for both 
removing (adding) electrons in symmetric states {k — 
0) and antisymmetric states c]^(fc = tt). 

In the strongly anisotropic case 1^3/2! = 1*1 3> \t'\ — 
1*5/2! = I ^1/2 1 1 the wave function of the Kramers' degen- 
erate ground state with Jz = ±15/2 has the simple form 
(for = +15/2) 

1^)15/2 = ^ (4/2(«) + ^4/2(^)) X 



^5/2 



{a)c\,^{a)cl,^{b)c\,^{h) 



(22) 



describing an independent jz = 3/2-electron in a mag- 
netically polarized background. This state is separately 
an eigenstate of the kinetic energy as well as the Coulomb 
energy. Removing an electron in the single-particle state 

"71 (''3/2('^) ~^ jf| 4/2 (^)) adding one in its or- 

thogonal counterpart leads to final states J(a) = A — 
J{h) and J(a) = 9/2 = J(b), respectively, where both 
sites are in the ground state multiplets of the local p and 
P configurations. As a consequence, we find two sharp 
peaks with full spectral weight at the energies oj = ±\t\ 
which evolve into a quasiparticle band in extended sys- 
tems. The additional spectral weight appearing in Figure 
Hat positive energies lo = +\t\ accounts for the prob- 
ability of adding an electron in the minority channel 
jz — —3/2. This generates final states which contain 
contributions from excited states like unfavorable local 
/^-configurations and excited multiplets in the local Z'^- 
subspaces. The overlap with states from the local ground 
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state manifold gives rise to a peak at w = +\t\ whose 
reduced height, however, highhghts the transfer of signif- 
icant spectral weight to the high-energy regime. Finally, 
the absence of dispersion in the jz = 1/2 and jz = 5/2 
spectral functions as well as the reduced spectral weight 
in the low-energy regime reveal the localized character of 
the corresponding orbitals in the ground state. The evo- 
lution of the spectra and in particular the reduction of 
spectral weight in the low-energy regime is reminiscent 
of the gap formation at the Mott-Hubbard transition. 
The qualitative difference between the 5/2- and the 1/2- 
channel as well as the asymmetry within the 1/2-channel 
reflect the intra-atomic (Hund's rule type) correlations. 

The structure of the spectral functions in the comple- 
mentary anisotropic limit 1^5/2! = l^i/z] — I'^'l ^ l^'l — 
1*3/2! can be explained by similar considerations. The 
Kramers' degenerate ground state is built from conflgu- 
rations where all but the jz — —3/2 {jz = 3/2) orbitals 
are occupied. This implies = ±3/2. It is important 
to note that this implies antiferromagnetic correlations 
between the sites. In this state, delocalization of the 
jz — 3/2 (jz = —3/2) is not principally excluded. It 
is, however, associated with breaking Hund's rules. As 
a result we find dispersion in the corresponiding chan- 
nels. The spectral weight, however, is reduced due to 
intra-atomic excitations. 

The characteristic structure of the spectra as summa- 
rized above is present also in the three-site cluster as can 
be seen from Figure The coherence effects, however, 
are more pronounced. For a discussion of the (approxi- 
mate) ground states in limiting cases we refer to )12 |. 

A striking feature is the splitting of the lower peaks 
for jz = ±1/2 and jz = ±5/2 at ii/2 = 1^/2 = -0.09, 
*3/2 ~ —0.01. The key to the explanation is provided by 
the observation that the corresponding ground state is 
characterized by antiferromagnetic intersite correlations, 
a property encountered already in the two-site cluster 
discussed above. The implications for the three-site sys- 
tem, however, are more subtle. For periodic boundary 
conditions the ground stae is a complex superposition of 
several (almost) equivalent states refiecting frustration. 
Removing an electron from the ground state leads to 
rather complex final states which, inturn , give rise to 
the observed splitting. 

The most prominent feature of the spectral functions 
in the strongly anisotropic limit is the partial localiza- 
tion which can be viewed as an orbital-selective Mott 
transition. The intra-atomic correlations strongly en- 
hance the anisotropics of the kinetic energy as can be 
seen from Figure^l The calculated projected kinetic en- 
ergy Tj./tj^ (Figure 21 ) agrees well with results from 
ground state wavefunction 12]. The small differences to 
the results from the ground state wave function are due 
to Lorentzian broadening with 77 = 0.005. 




0.05 0.1 0.05 0.1 



Figure 5: Phase diagrams of the three-site cluster filled with 
8-electrons for negative hopping parameters tj^ , derived for 
the total magnetization Jz for periodic boundary condition 
(left panel) and open boundary conditions (right panel). 



V. CLUSTER PERTURBATION THEORY 

The essential benefit gained from CPT is the possibil- 
ity to calculate the spectral functions at arbitrary wave 
vectors and band filling. In the present paper we explore 
the first aspect treating CPT as a sophisticated inter- 
polation scheme for the variation with k. The averaged 
particle number per site, on the other hand, is kept fixed 
at the value used in the underlying cluster calculation. 
Comparing the spectra of the infinite chain to those of 
isoelectronic finite clusters provides us with an estimate 
of finite size effects since the building blocks of the chain 
are treated with different boundary conditions. 

The qualitative behavior of the spectral functions re- 
fiects the detailed structure of ground state which is de- 
termined by the competition between the energy gain 
due to intra-atomic correlations and kinetic energy. The 
latter, however, depends upon the boundary conditions 
imposed upon the wave functions. The cluster calcula- 
tions discussed in the preceding section were performed 
adopting the usual periodic boundary conditions. The 
CPT, however, requires cluster Green's function for open 
boundary conditions. In the strong-coupHng regime, the 
dependence upon the transfer integrals of the ground 
state is only weakly affected by the actual choice of the 
boundary conditions. In the strongly anisotropic limit 
with one dominating transfer integral (a) the high-spin 
states with ferromagnetic inter-site correlations are ener- 
getically most favorable. They are followed by rather 
complex intermediate spin states (b) as the case of 
isotropic hopping (c) is approached. Along this line the 
salient feature of the correlated ground state is its high 
degeneracy. In the case of two dominating hopping chan- 
nels (d) low-spin phases with anti-ferromagnetic inter-site 
correlations are formed. These features originally derived 
for periodic boundary conditions are encountered also for 
open boundary conditions as can be deduced from Figure 

m 

Figure displays the orbital-projected spectral func- 
tions of the linear chain. The data refer to the low- 
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Figure 6: Spectral functions of a one-dimensional 5/-system 
from the CPT with cluster of two sites (upper panel) and 
three sites (lower panel). The orbital-projected spectra are a 
superposition of the +jz and ~jz parts of the spectra. The 
anisotropic hopping parameters are and t3/2 = —0.09 eV and 
= ^5/2 = —0.01 eV . The peaks are broadened with a 
finite imaginary part of = 0.03. 



Figure 7: Spectral functions of a one-dimensional 5/-system 
from the CPT with cluster of two sites (upper panel) and 
three sites (lower panel). The orbital-projected spectra are a 
superposition of the +jz and —jz parts of the spectra. The 
anisotropic hopping parameters are and = —0.01 eV and 
ti/2 = ^5/2 = —0.09 eV . The peaks are broadened with a 
finite imaginary part of = 0.03. 



energy regime where the electronic excitations involve 
transitions between the lowest Hund's rule multiplets. 
We adopt the same values for the transfer integrals as 
in the cluster calculations. The system is clearly in 
the strong-coupling regime where the symmetries of the 
ground states are rather sensitive to the detailed descrip- 
tion of the kinetic energy. The high-energy features as- 
sociated with transitions to excited atomic states are not 
shown here. 

As regards the variation with wave vector of the spec- 
tral functions, the qualitative behavior closely parallels 
the one obtained for the finite clusters. The spectral 
functions clearly show a dispersing band in the orbital 
channel with the dominant hopping. The position of the 
narrow peak whose finite width is to be attributed to the 
additional broadening follows the characteristic cosine 
dispersion of a one-dimensional tight-binding band with 
nearest-neighbor hopping. The gap in the dispersion of 
the dominant hopping channel in Figuredare introduced 
by the antiferromagnetic intersite correlations which ef- 
fectively reduce the size of the Brillouin zone. The band 
width, however, is reduced by a factor of ~ 7/12 as com- 
pared to 4 \t\ = 0.36eV expected for uncorrelated elec- 
trons. The renormalization reflects the transfer of spec- 
tral weight to the high-energy satellites. The line widths 
of the dispersing bands in the linear chain are approxi- 
mately twice the value of the bonding-anti-bonding split- 
ting in the two-site cluster. This discrepancy is due to 
the fact that the number of nearest neighbors in the lat- 



tice is twice that in the corresponding cluster. The 5/ 
channels with subdominant hopping exhibit low-energy 
excitations which are incoherent and which may exhibit 
a pseudo-gap at the Fermi energy. It should be noted 
that the CPT spectra based on three-site clusters are in 
good agreement with their two-site cluster counterparts. 
The cluster size and the boundary conditions apparently 
affect the incoherent excitations in the channels with the 
subdominant hopping. This can be seen from compar- 
ing the CPT spectra obtained from two- and three-site 
clusters. 

The co-existence of a coherent 5/ band and incoherent 
/-derived low-energy excitations implies partial localiza- 
tion. The orbital-selective localization suggests the pres- 
ence of different types of low-energy excitations which 
have different orbital character. Photoemission experi- 
ments with polarized incident light should be able to dis- 
tinguish the two. In fact, recent experiments on URu2Si2 
seem to be consistent with this hypothesis. 

The Momentum Distribution FXmction (MDF) as de- 
fined in Eq. © is displayed in Figure |HI The variation 
with wave vector k clearly shows the dual character of 
the 5/ electrons. In the homogeneous high-spin phase, 
the height of the discontinuity in Uj^ (fc) for orbitals with 
dominant hopping is consistent with the spectral weight 
of the dispersive low-energy peak. In the low-spin phase, 
however, it reflects the reduction of the Brillouin zone 
due to antiferromagnetic correlations. 

Figure El displays the orbital-projected DOS as calcu- 
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Figure 8: Orbital-projected MDF calculated within CPT 
starting from three-site cluster and the hopping parameters 
^3/2 = t and ti/2 = = t' . The spectral functions displayed 
in Figure were integrated over energy in an interval of —4 eV 
below the chemical potential. 
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Figure 9: Orbital projected density of states Dj^ (a;) of a 
one-dimensional 5/-system from the CPT with cluster of 
three sites. The anisotropic hopping parameters are and 
t3/2 = -0.09 eV and ti/2 = ^5/2 = -0.01 cV . The peaks 
are broadened with a finite imaginary part of q = 0.03. 



lated within CPT. 



trons by applying CPT. The result derived for a one- 
dimensional chain confirm the idea of orbital-dependent 
localization in anisotropic systems. The orbital-projected 
spectral functions for the channels with dominant hop- 
ping display narrow dispersive bands in the vicinity of 
the Fermi energy while the low-energy single-particle ex- 
citations in the remaining channels involve mainly in- 
coherent local transitions. In all channels, considerable 
spectral weight is transferred to the high-energy regime 
where we find intra-atomic transitions to excited mul- 
tiplets as weh as valence transitions. The transfer of 
spectral weight rcnormalizcs the widths of the dispers- 
ing bands and hence increases the corresponding effective 
masses. 

The dual character of the 5/ electrons, i. e., the di- 
vision of 5/ states into delocahzed and locahzed orbitals 
can be clearly seen form the orbital-projected momentum 
distribution function. The latter has a pronounced dis- 
continuity for the orbitals with dominant hopping while 
its smooth variation can almost be neglected in the re- 
maining channels. In a three-dimensional crystal such 
a behavior implies that only the orbitals with dominant 
hopping contribute to the quantum oscihations of the 
dHvA effect. 

Calculations starting from clusters with two and three 
sites yield similar results in the strong coupling regime. 
The orbital-dependent suppression of the kinetic energy 
calculated for the chain agrees well with previous re- 
sults obtained for small finite clusters imposing periodic 
boundary conditions. These finding suggest that the 
features associated with orbital-selective locahzation are 
rather robust. 

The model calculations presented suggest that the 
orbital-dependence of the 5/ spectral functions can be 
used as fingerprint for orbital-selective localization in ac- 
tinide systems. The orbital dependence could be observ- 
able by varying the polarization of the incident light. 



VI. SUMMARY AND OUTLOOK 

We have studied the consequences of strong intra- 
atomic correlations on the spectral functions of 5/ elec- 



It is therefore desirable to perform similar calculations 
for models describing real materials. This requires the 
calculations to be extended to higher dimensions. An 
important generahzation is the inclusion of non-diagonal 
transfer integrals which is currently in progress. 
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